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ABSTRACT 

The first complete submillimetre spectrum (190-670/im) of the Seyfert 2 galaxy NGC 1068 has been 
observed with the SPIRE Fourier Transform Spectrometer onboard the Herschel Space Observatory. 
The sequence of CO lines (J up =4-13), lines from H 2 0, the fundamental rotational transition of HF, 
two o-H 2 + lines and one line each from CH + and OH + have been detected, together with the two 
[CI] lines and the [NII]205/xm line. The observations in both single pointing mode with sparse im- 
age sampling and in mapping mode with full image sampling allow us to disentangle two molecular 
emission components, one due to the compact circum-nuclear disk (CND) and one from the extended 
region encompassing the star forming ring (SF-ring). Radiative transfer models show that the two 
CO components are characterized by density of n(H2)=10 4 ' 5 and 10 2 9 cm~ 3 and temperature of T^ n 
=100K and 127K, respectively. The comparison of the CO line intensities with photodissociation 
region (PDR) and X-ray dominated region (XDR) models, together with other observational con- 
straints, such as the observed CO surface brightness and the radiation field, indicate that the best 
explanation for the CO excitation of the CND is an XDR with density of n(H 2 ) ~ 10 4 cm -3 and 
X-ray flux of 9 erg s _1 cm~ 2 , consistent with illumination by the active galactic nucleus, while the 
CO lines in the SF-ring are better modeled by a PDR. The detected water transitions, together with 
those observed with the Herschel PACS Spectrometer, can be modeled by an LVG model with low 



temperature (Tki n ~ 40K) and high density (n(Ha) in the range 10 



,6.7 



10 cm 3 ). The emission of 



H2O" 1 " and OH + are in agreement with PDR models with cosmic ray ionization. The diffuse ionized 
atomic component observed through the [NII]205/im line is consistent with previous photoionization 
models of the starburst. 

Subject headings: Galaxies: individual: NGC 1068 - Galaxies: ISM, nuclei, active, starburst, Seyfert 
- Techniques: imaging spectroscopy 



1. INTRODUCTION 



NGC 1068 (Messier 77) is a nearby (cz = 1 137 km 
1 ^ and bright (Lfp 



-1000/^m 



-2 x 10 11 L, 



O: 



Bland- 



Hawthorn et al.|1997[) Seyfert galaxy, often considered as 
the prototypical Seyfert type 2 galaxy. However, since 
the discovery of the broad permitte d lines in the polar- 
ized optical spectrum of NGC 1068 ( |Antonucci k, Miller 
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1985 1, it has become clear that this galaxy was in reality 
a hidden broad line region galaxy. The general distinc- 
tion between the two types of Seyfert galaxies might be 
due only to orienta tion effects, acco rding to the so called 
Unification model ( Antonucci]|1993[ ) . 

Being the strongest nearby Seyfert 2 galaxy, it has 
been observed extensively over the whole electromagnetic 
spectrum. Molecular (CO and HCN) observations have 
shown a prominent starburst ring (hereafter SF-ring) at 
a radius of 1.0-1.5 kpc and a central circum-nuc lear disk 
(hereafter CND) with a diameter of d ~ 300 pc ( Tacconi 
et al.|1 994; Sc hinnerer et al.|2000|). Near-IR observations 
( |Scoville et al.|[l988[ |Thronson et al. 1989) have clearly 
revealed a 2.3 kpc stellar bar. The compact (~1 pc) hot 
dust source in the nucleus of NGC 1068, measured us- 
ing near-i nfrared speckle ima ging and integral field spec- 
troscopy ( Thatte et al.|[l997 ), is probably heated by the 
AGN's strong radiation field and possibly associated with 
the postulated dense circum-nuclear torus. A jet was 
observed from centimeter to millimeter wavelengths ex- 
tending out to several kiloparsecs from the center (e.g., 
Krips et al.||2006| |Gallimore et al.||2004| . Mid-IR obser- 



vations revealed hot and ioniz ed~g as bico nically following 
the path of the radio jet (e.g., Miiller Sanchez et al.|2009 



and references therein) and inTJIcating the existence of a 
parsec-scale warm dust torus ( Jaffe et aL||2004 ). Recent 
interferometric observations 01 (JO (3-2) and CO (1-0) 
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show that the CO (3-2) emission peaks in the central re- 
gion within ~ 5 " from the nucleus, while the CO (l-O) 
emission is mainly located along the spiral arms (Tsai 
eTaI][2012| . 



Spectroscopic coverage of NGC 1068 has been exten 
sive. In the mid-IR to far-IR wavelength range, spectr; 



tiave been measured by the ISO ( Kessler et al. 1996 


) SWS 


(de Graauw et al. 


19961) and LAVS (Clcgg cfani 


|l996fc 


spectrometers (Lutz et al. 


2000 
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covering 

spectral ranges. At millimeter wavelengths it was re- 
cently observed from the ground in the 190-30 7 GHz 



(976-1578 /im) range (Kamenetzky et al. 2011). The 



submillimctcr waveband is one ot the tew spectral re 
gions that have not been so far explored; the new obser- 
vations made by t he Spectral and Ph otometric Imaging 
Receiver (SPIRE) ([Griffin et afl2 010) Fourier Transform 
Spectrometer (FTS) ( |lNaylor et al. 2010b ) onbo ard the 
Herschel Space Observatory ( Pilbratt et al.|2010 |, cover- 
ing the spectral range from 190 /im to 670 /im, till most of 
this gap with the first complete submillimeter spectrum 
of NGC 1068. 

Submillimeter spectral measurements are of particular 
interest as NGC 1068 is a prime candidate to study the 
effects of the AGN onto the circum-nuclear material and 
the surrounding disk. In particular, the study of the ex- 
citation conditions and chemistry of the CND appear al- 
ready from the existing ground-based molecular line ob- 
servations to be very different from the starburst galaxies 
environments. The peculiar line ratios of different molec- 
ular transitions, mostly HCN, HCO+, and 12 CO, led to 
the suggestion that the CND of NGC 10 68 harbors a gi- 
ant X-ray-dominated region (XDR, e.g., fRotaciuc et al. 
19911 |Maloney et al.|1996l|Usero et al.|2004||Kohno et al. 



20081). The HCN and HC O + molecular line studies of 
Krips et al. (2008 2011) confirm an increased abun- 
dance ot HCN and/ or increased kinetic temperatures. 
CO lines can also discriminate between "classical" pho- 
todissociation regi ons (PDRs) and X-ray do minated re- 
gions (XDRs) (e.g. |Meijerink fc Spaans|2005|). Using the 
intermediate J rotational lines from the CO molecule, 
from J up =4 to J up =13, which can be observed with the 
SPIRE FTS, we want to test if indeed in NGC 1068 an 
XDR is needed to explain the spectral line energy distri- 
bution originating from the CND. The case of the ultra- 
luminous IR galaxy Mrk231 already demonstrated that 
the Herschel-SPIRE data are indeed able to discrimi- 
nate between the two emission me chanisms and there- 
fore detect the effects of the AGN (van der Werf et al. 
Herschel-PACS observations ot the high-J CO 
> 14) detected in NGC 1068 have been re- 



20101. 
lines (J up 

cently presented by |Hailey-Dunsheath et al. (2012). The 
two components, at high and medium excitation, needed 
to explain the observed CO lines arising from the central 
10" region from J up = 14 to J„ p = 24 can be excited by 
X-ray or shock heating, while far-UV heating is unlikely. 

The _ffer.se/ie/-SPIRE spectroscopic observations of 
NGC 1068 presented here have been collected under the 
guaranteed time key project "Physical Processes in the 
Interstellar Medium of Very Nearby Galaxies" (PI: Chris- 
tine Wilson). Within the same observational program, 
SPIRE and Ph otodetector Array Ca mera and Spectrom- 



eter (PACS) (]Poglitsch et al. 2010[ ) photometric images 
have been collected, which will be presented in a forth- 



coming paper with the detailed analysis of the continuum 
emission (Spinoglio et al 2012, in prep.). 

2. OBSERVATIONS 

2.1. SPIRE Spectroscopy 
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Fig. 1. — Observed source sizes, as measured from the mapping 
observations of the brightest SPIRE lines (FWHM in arcsec, with 
1<T error bars) as a function of frequency (colored points) compared 
with the SPIRE-FTS beam (broken lines). All source sizes in the 
SSW, except for [Nil] 205 (im, are consistent with the beam width 
at that frequency. In contrast, the source sizes in the SLW are 
more extended than the beam. 



NGC 106 8 was observed with the SPIRE-FTS .Griffin 
et al.|2010 ) onboard the Herschel Space Observatory ( |Pil" 
bratt et al.| |2010 ) both in the single pointing mode with 



sparse image sampling and in the mapping mode with 
full image sampling. The FTS has two detectors arrays 
called the spectrometer long wave (SLW, in the range of 
303-671 /im) and the spectrometer short wave (SSW, in 
the range 194-313 /im), with a small (10 /im) overlap in 
wavelength. The spectral resolution of the FTS ranges 
from about 280 km s _1 to 950 km s _1 in the high reso- 
lution mode, moving in wavelength from 194 /im to 671 
/im. The two SPIRE FTS observations of NGC 1068 were 
collected on Operational Day 626 (29 January 2011); the 
fully sampled map has Observation ID 1342213444 and 
the single pointing "deep" spectrum has Observation ID 
1342213445. The pointed observations were carried out 
both at low and high spectral resolutions (FWHM ~ 30 
GHz and 1.44 GHz), covered with 32 repetitions each, 
for a total on-source inte grat ion time of 410 seconds and 
4,262 seconds respectivel}j 10 | The shallower mapping ob- 
servations at high resolution had 8 repetitions and cov- 
ered an area with a diameter of approximately 2'. The 
total on-source integration time was 17,050 seconds. 

The FTS beam size and shape vary as a function of fre- 
quency and cannot be characterized by a simple Gaussian 

10 See the Herschel SP IRE Observers manual, available at 
http: / /herschel. esac.esa.int/Docs/SPIRE/html/spire_om. html 
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Fig. 2.— Left: CO(3-2) integrated intensity map o f the centra l regio n of NGC 1068 from |Tsai et a,I.| (|2012[ >. 
intensity map of the central region of NGC 1068 fromlTsai et al.l (12012b . The minimum SSVV (blue) and the m 
FTS beam sizes of 17" and 42", respectively, are shown m the figures as concentric circles. 



Right: CO(1-0) integrated 



response. As can be seen from Figure [TJ where the ob- 
served source sizes of the brightest lines detected by the 
mapping observations of the SPIRE FTS are plotted as a 
func tion o f frequency (see the method described in Sec- 
tion [2TT|), the FTS beam widths (FWHM) range from 
17" to 42 , probing a wider area in the sky with increas- 
ing wavelength, and show a strong discontinuity between 
the SLW and SSW bands, due to the multi-moded feed- 
horns used for the spectrometer arrays (see Section 4.2.3 
and Figure 5.13 in the SPIRE Observers' Manua£°}. 

The combination of the complex morphology of the 
inner region of NGC 1068 with the SPIRE-FTS beam 
properties yield a particular coupling effect between the 
two. Figure [2] shows the CO (3-2 ) and CO(1-0) inter- 
ferometric maps from Tsai et al. ( |2012 ) of the central 
region of NGC 1068. As can be seen from these maps, 
the morphology of NGC 106 8 is dominated by two main 
components (see also, e.g., |Garcia-Burillo et al. |2010 
and references therein): the compact (of the order o 
4" in diameter) CND at the center of the galaxy, possibly 
associated with the active nucleus, and a relatively ex- 
tended ring, with a radius of the order of 10-20", whose 
emission is dominated by star formation (SF) activity. 
The relative contribution from these two spatial com- 
ponents to the various rotational transitions of the 12 CO 
molecule strongly varies with the transition from the low- 
J (Jup ~ 4) to the high-J (J up ~ 13) lines. In this range 
the dominant source of CO excitation changes from the 
extended SF-ring to the compact CND. There is also a 
variation of the angular scale probed by the FTS beam, 
from extended (beam FWHM ~ 42") to compact (17"), 
as shown by the two concentric circles in the maps of 
Figure [2} 

2.1.1. Source size measurements 

We first describe the reduction procedure tha t was u sed 
for so urce size measurements, while in Sections |2 . 1 . 2| and 
|2.1.3| we describe the reduction procedure for spectral 
line measurements. We employ two reductions as the 
two measurements have different requirements. 

The mapping observations of NGC 1068 have been re- 
duced with the standard reprocessing pipeline in HIPE 
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Fig. 3. — The observed submillimeter spectrum of NGC 1068, as 
measured from the SPIRE-FTS pointed observations. As can be 
seen the two SSW and SLW spectra do not overlap. This is due to 
both the different beam sizes and the intrinsically different source 



v.7 and converted into spectral cubes using the Naive- 
Projection tasl^with custom pixel sizes of 11" and 
13.8" for the SSW and SLW bands, respectively. These 
sizes were chosen empirically starting from a large size 
and reducing it until holes (pixel with zero coverage) ap- 
peared in the maps. The final cube for each of the two 
FTS bands has a variable on-source integration time in 
the map, between 16 and 32 FTS scans. From these 
two cubes a number of sub-cubes have been obtained 
by extracting a 60 GHz wide window (corresponding to 
~ 200 bins in high spectral resolution) around each line 
detected in the spectrum of the brightest spaxel, where 
a spaxel is defined as each one of the 5x5 spatial pixels 
which fill the spectrometer field of vie , wp 2 ] and modeling 
and subtracting the continuum emission with a I s * order 



11 This and 
scribed in the 
available at: 



the following routines are de- 
SPIRE Data Redu ction Guide, 
http: / /herschel.esac.esa.int/hcssdoc- 
6.0/load/spire_dum/html/spire_dum.html 

See the FA(JS Ubserver's Manual at: 
http: / /herschel.esac.esa.int/Docs/PACS/html/pacs_om.html 
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Fig. 4. — The observed submillimetre spectrum of NGC 1068, as measured from the SPIRE-FTS pointed observations. From top to bottom 
are presented the consecutive spectral ranges. The emission line positions are indicated. Line identifications are given for detections with 
S/N>3tr (see Table [l}. 



polynomial using HIPE's task LineIntensityMap. 

The FTS data from a pointed observation have two 
available calibrations which provide the target's correct 
flux for the two extreme cases of an unresolved source and 
of an extended source that completely and uniformly fills 
the instrument beam. To determine which calibration to 
use, we measured the emitting region size for several lines 
from the analysis of the spectra cubes (see Figure [l}. 

The source sizes of all the transitions present m the 
SSW band, with the exception of the [Nil] 205 /jm line, 
are consistent with originating from an unresolved source 
and they are therefore well calibrated in the pointed spec- 



trum using the unresolved source solution. For this rea- 
son, we used the line fluxes measured from the pointed 
spectrum for all the transitions in the SSW band except 
for the [Nil] 205 /zm line, for which we have used the 
mapping observations. ([Nil] emission arises from ion- 
ized gas and so we may expect some contribution from 
the SF-ring.) For all the transitions in the SLW band, 
we measured the source sizes from the mapping observa- 
tions. 

We have obtained line integrated maps with the 
task IntegrateMapFromCube from the spectral sub- 
cubes. The integration over the lines has made no use of 
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the fitting functionality inside the task, because of the 
possible spurious detections observed for spectral fea- 
tures with low signal-to-noise ratios. These line maps 
have been used to measure the angular extent of the re- 
gion emitting every detected line using a 2D Gaussian 
fit. While the instrument PSF is not strictly gaussian, 
because a large fraction of its light is in wider secondary 
lobes, this approach is robust enough to quantify the 
scale size of the emitting region through the fit of the 
PSF core. 

In conclusion, the deep pointed spectrum in the SSW 
band samples primarily, not exclusively, the lines emit- 
ted by the compact source, while the SLW band mea- 
sures lines emitted by both the compact and the extended 
structures. This coupling effect between the beam and 
the target morphology makes it difficult to reconstruct 
the correct spectral line energy distribution (SLED) for 
the 12 CO molecular lines for the CND alone. To over- 
come these problems, we have combined the information 
extracted from the pointed spectrum with those derived 
from the mapping observations. 

2.1.2. Spectroscopic mapping observations 

The mapping observations of NGC 1068 were again re- 
duced using a modified version of the standard Spectrom- 
eter Map pipeline script in the Herschel Interactive Pro- 
cessing Environment (HIPE) v. 9, developers build 588 
and version 8.1 of the SPIRE calibration context. The 
standard reduction assumes that the source is extended 
and uniformly fills the beam. Since our source does not 
fill the beam, the unique point source flux conversion 
measured for each bolometers was applied to each of the 
bolometers in our two arrays. After the point source 
correction, spectral cubes for the SSW and SLW were 
created using the spireProjection task with projection- 
Type— "naive" with a custom pixel size of 11" and 13.8" 
respectively as before. This cube was used for all SLW 
and the [NII]205/im line intensity measurements. 

2.1.3. Pointed observations 

The single pointing/sparse observations of NGC 1068 
have been reduced using the Spectrometer Single Point- 
ing user pipeline in the Herschel Interactive Processing 
Environment (HIPE) v. 9, developers build 588 and ver- 
sion 8.1 of the SPIRE calibration context. An older ver- 



sion of this pipeline is described in |Fulton et al.| ( 2010[ ). 
The observed submillimeter spectrum of 1NGC1068, as 
measured from the SPIRE-FTS pointed observations, is 
shown in Figure [3j where the large offset at ~ 300/iin is 
due to the large difference in the beam of the SSW and 
the SLW spectrometers. 

2.2. Line intensities 

As explained in the previous section, the transitions 
measured in the SSW band, except for the [Nil] 205 fim 
line, have been considered spatially unresolved and there- 
fore measured from the pointed/deep spectrum, while, on 
the contrary, the other partially extended lines have been 
measured from the mapping observations. 

In order to measure the line fluxes of the extended 
compon ent, the line maps have been convolved with 2D 
kernels (Bendo et al. 112012) to reach the same beam size 
TEe~fi: 



coarsest angular resolution (FWHM = 42") reached by 
the FTS for a 12 CO line (i.e. the CO (4-3) at 461.04 
GHz). The kernels were built using the SPIRE-FTS PSF 
measurements, dividing the PSF image of the CO (4-3) 
line by that of each detected line in the Fourier space. 

The baseline was removed by first masking all of the 
spectral lines before fitting a high order polynomial to 
the remaining continuum. The intensities of the detected 
lines have been computed with a Levenberg Marquardt 
fitting procedure using the sine function model. Table [T] 
presents, for each detected transition, the rest frequency, 
the energy of the upper level (in K) , the beam FWHM (in 
arcsec.)(see the SPIRE Observers' Manuap°| the inten- 
sity value with both the la statistical uncertainty and 
total uncertainty, which include the calibration uncer- 
tainty added in quadrature. For the mapping observa- 
tions, we have adopted a fixed FWHM of 43.4" for all 
detected lines. For the pointed observations of the com- 
pact source, we have adopted a calibration uncertainty 
of 10%, while for the mapping observations of the ex- 
tended source, in addition to the calibration uncertainty, 
we have also added in quadrature a flat-fielding uncer- 
tainty of 7% (see the SPIRE Observers' ManuaQ- Th e 
lines detected with the FTS are not spectrally resolved. 

Figure [4] presents the observed spectrum of NGC 1068 
with the identifications of all lines detected at S/N>3. 
A total of 24 emission emission lines have been detected 
with the Herschel-SPIRE FTS: ten lines from CO (from 
Jup=4 to J up =13), six lines from H 2 0, two from H 2 + , 
one from OH + , HF and CH + each, in addition to the 
three atomic lines (from [Nil] and [CI]). The apparent 
ringing in the spectrum is due to the sine response of 
the FTS to each line. 

In Table m we also included literature measurements, 
relative to both the compact source and the extended 
source. In particular, we considered in our analysis the 
interfero metric observations of CO (1-0), (2-1) and (3- 
2) from Krips et al. (2011), and the water lines from 
the Herschel-PA<Jti observations (S. Hailey-Dunsheath 
2012, private communication) for the compact compo- 
nent, and, for the extended component, CO (2-1) Z-Spec 
spectrometer observ ations at the Caltech Su bmillimeter 
Observatory (CSO) ( Kamenetzky et al.|2011 ) and James 
Clerk Maxwell Telescope (JUIVIT) (JU (3-2) observations 
convolved to the 43.4"beam (C. Wilson 2012, private 
comm.). 

3. RESULTS AND DISCUSSION 

3.1. CO Radiation transfer modeling 

To obtain a general solution for the CO physical con- 
ditions, we have first used the radiative transfer code 



at every frequency. 



final PSF corresponds to the 



RADEX ( vanderTaket al.|2007[ ). Then, considering the 
best possible ranges of the parameters as derived from 
RAD EX, we have also u sed the PDR and XDR models 
fro mlWolfire et al.|(|2010| a nd |Meijerink fc Spaansj ( |2005 1 
and Meijerink et al. "(2007), respectively, to provide con- 
straints on the origin of the excitation of the measured 
lines, and in particular to determine whether there is an 
influence of the AGN, through its X-ray emission, on the 
molecular gas in the CND and SF-ring associated with 
NGC 1068. 

3.1.1. RADEX models 
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RADEX is a non-Local Thermodynamic Equilibrium 
(LTE) code available fro m the Leiden Atomic and Molec- 
ular Database (LAMDA, |Schoier et al.|2005[ ). Under the 
assumption of a uniform medium, using an escape prob- 
ability formalism that models the entire emitting region 
as a single zone, RADEX performs statistical equilib- 
rium calculations involving collisional and radiative pro- 
ce sses. We use d the c ollisional rate coefficients with H2 
of Yang et al. (2010). Considering the radiation field 
from background sources given as input, it computes the 
molecular level populations in the optically thin limit in 
the first iteration and the optical depths and the escape 
probabilities for each line. The interdependence of the 
molecular level populations and the local radiation field 
requires that the solution is reached through an iterative 
method. At each iteration, the code computes new level 
populations using new optical depth values, until the two 
converge on a consistent solution. Then the resulting line 
intensities are given as output. 

The input parameters required for RADEX are the 
background radiation field, the kinetic temperature Tki n , 
the molecular hydrogen number density n(H2) (assumed 
to be the only collision partner), and the ratio between 
the column density of the molecule and the width of the 
lines, which is the physically relevant quantity that de- 
termines the optical depth. 

We used RADEX to create a grid of solutions for the 
CO transitions spanning wide ranges in T^m, n(H2), and 
Nco ( see Table [2]) . To compare this grid of models with 
our CO meas urements we used the Bayesian likelihood 
analys is code (|Ward et al.|2003[), developed by t he Z-Spec 
Team qjNayfor et al.||2010a[ |Kamenetzky et al.||2011[ ). 
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Fig. 5. — RADEX model fitting for the compact component as- 
sociated with the CND. Square points (red) represent the PACS 
data, circle points (blue) depict the FTS deep spectrum data, and 
the diamo nd points (grey) sh ow the interferometric observations of 
the CND (Krips et al. 2011), the dotted line shows the ME model 
of |Hailey-D"u nshcath et al. (2012), while the dashed line shows the 
best fitting configuration 01 KADEX parameters (see column 4D 
Max of Table [3). The sum of the two models is shown as a solid 
line. 



Computing the RADEX models, we have considered 
both models with and without the inclusion of the lo- 
cal background radiation to be added to the cosmo- 
logical component of the CMB at 2.73K. To estimate 
the local background in the far-IR, we have fitted the 
aperture photometry of the Herschel-PACS and SPIRE 
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Fig. 6. — RADEX model fitting for the extended component as- 
sociated with the SF-ring. Notations as in Figure [5] The crosses 
(yellow) represent the FTS measurements from the spectral map 
and the white stars show the ground-based o bservations of the ex- 
tended component |[Kamcnctzky lit al.|201l| Wilson 2012, private 
communication). The dot-dashed line shows the best fitting con- 
figuration of RADEX parameters (see column 4D Max of Table[3[l. 
The thin solid line is the best model for the CND component. Tne 
sum of all models is shown as a thick solid line. 



maps at 70, 160, 250, 350 and 500 fim, respectively, that 
will be presented and discussed in a forthcoming paper 
(Spinoglio et al. 2012, in prep.). We have used a circular 
aperture with a diameter of 33" and a dust emissivity 
law with /3=2 and we have derived a gray-body temper- 
ature of 35 K. In agreement with the h i gher J -CO line 
fitting results of Hailey-Dunsheath et al. (2012), we have 
found that the inclusion of this local background does 
not change significantly the results of the CO fits that 
only include the CMB. Therefore we have adopted only 
the CMB black-body radiation at 2.73 K as background 
radiation, for both the extended and compact sources. 

The presence of several spatial and spectral compo- 
nents in the CO SLED has required a multi-step ap- 
proach, looking for the best solution of one component at 
a time and subtracting it from the others. We started by 
considering the CO SLED of the compact source, com- 
posed of the FTS-SSW measurements in combination 
with the compact interferometric observations. For the 
high-J (J up > 14) PACS data, we used the medi um exci 
tation (ME) component LVG model presented by Hailey 
Dunsheath et al. ( |2012 ), characterized by a density of 



10 a Y cm _J and a kinetic temperature T^in = 150 K. The 
ME model for the J up > 13 CO lines was subtracted from 
the observed CND spectrum before attempting to fit the 
SLED. 

In Table [3j the results obtained with the Z-Spec code 
are summarized. For each of the two models of the CND 
and SF-ring the following parameters are given: the ki- 
netic temperature of the CO, Tk; n , the molecular hydro- 
gen density, n(H 2 ), the CO column density, Nco, the 
area filling factor, $A: the gas pressure, P, the beam 
averaged CO column density, (Nco)j and the molecular 
hydrogen mass M(H2), that has been computed using a 
CO a bundance of xgo = 3 x 10~ 4 , adopted by the Z -Spec 
code ( Naylor et al.| 2010a Kamenetzky et aL]|2011 1. 

For each parameter, the median value, the l-a interval, 
the ID and the 4D maximum values are given. ID refers 
to the maximum value of the integrated parameter dis- 
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tribution while 4D refers to the value of that parameter 
at the best fit solution. 

The compact source is unresolved with all the consid- 
ered instruments and identified with the CND, while the 
extended source data probe the emission arising from a 
region of 43.4" in diameter, compatible with the SF-ring 
structure. We have obtained the best fit to the CND with 
RAD EX, as can be seen in Figure [5j by fitting the mid- J 
lines (from J„ p = 9 to J up = 13) together with the low- 
J lines observed with interferometric techniques, which 
isolate the CND component, and subtracting the contri- 
bution from the medium excitation (ME) compo nent of 
the LVG model of |Hailey-Dunsheath etaL] ( |2012[ ). 

The best fit of the SF-ring with RADEX (seeTigure|6| 
was then obtained by subtracting from the observed lines 
(from J up = 4 to 3 up = 8 with SPIRE and 3 up = 3 and 2 
with JCMT) the contribution arising from the compact 
source by summing the mid- J and high- J contributions. 
As can be seen in Figure |6j the extended source gives 
a significant contribution to the over-all SLED only for 
J«p <9. To summarize, as can be seen in Table [31 we 
found that the best fit RADEX model results for the 
CND give T kin =100 K and n(H 2 )= 10 4 5 cm" 3 and for 
the SF-ring T kin =127 K and n(H 2 )= 10 29 cm- 3 . 

3.1.2. PDR and XDR models of the CND component 



14 to 



presented in Table [lj higher J lines (from J„ p 
J MD = 24) have been observed b y PACS and presented in 
Hailey-Dunsheath et al.| ( 2012 1. We have tried to model 
the whole sequence of CO lines detected with PDR/XDR 
models, however no single component was able to fit the 
data; three distinct components were needed, all origi- 
nating from the compact CND, to fit the CO lines down 
Moreover, as outlined in Hailey-Dunsheath 
in principle 



to Jut) — 9- 



et al.[ ( |2012 |), even if 
fit the data 



PDR models could 



'rom J up = 14 to J up = 24, the morphology of 
the H 2 near-IR emission and the two different kinematics 
of the CO lines at J„ p < 17 and 3 up > 20, indicate that 
these two sets of lines trace physically distinct compo- 
nents, therefore excluding a PDR origin of these high-J 
CO lines. Consequently, for the high- J (J uv > 14) PACS 
data we used the model presented by |Hailey-Dunsheath| 



et al. (2012), consisting of two XDR components with 



cm 

-2 



3 and incident X-ray fluxes of 9 
respectively. Then we have per- 



a density of 10 5 75 
and 160 ergs -1 cm 
formed the analysis on both the compact source and the 
extended source, once the other components were fitted 
and removed as explained in the previous section. 

The mid-J lines originating from the CND are best de- 
scribed by a sequence of PDR models ranging from an in- 
cident FUV flux of G - 10 5 and density of 10 45 cm -3 , 
in agreement with the RADEX model, to a FUV flux of 



PD R and XDR models f r om [W olfire ^Tl (pOTol) and £«, ~ ^ ° f ,~ T'^nlT ^ 

r ■■ — ~i — p — c 1 hinhzk, j i ii"' 1 '■ I ' , — n figure \7n. On the other hand, also XDR models either 

from Meijermk & Spaans (2005) and Mciicnnk et al. ° U m , ,„«n , • . , 

T^ i i u i i — rr — 1 i — j 1 ■— : — n , , t 1 with density ~ lO^cm or ~ 10 D U cm J and incident 



(2007P 3 describe the thermal and chemical balance of 
molecular gas that is exposed to far-ultraviolet (FUV) 
radiation (6-13.6 eV) and X-rays (1-100 keV), respec- 
tively. The codes used to determine the gas conditions 
in these regions as a function of depth take into con- 
sideration elaborate chemical networks and the cooling, 
heating and chemical processes which are completely de- 
termined by the radiation field. 

We have c onsidered the grids of PDR models of Wolfirc 
et al. (20101, that span from 10 1 to 10 7 cm -3 in density 
and from lO" 5 to 10 6 5 G tFI in FUV flux The PDR 
models are based on those of Kaufman et al.|fl2006|), with 



updat es from |Wolfire et al. ( |2010p and |Hollcnbach et al. 
(2012). In particular th e atomic and molecular fr eeze-out 



and grain chemistry in Hollenbach et al. (2012) are in- 
cluded, as well as their PAH photo rates which affect the 
ion- neutral chemistry and the production of CO. We use 
their "linear y ield" PAH° p hotoionization rate slightly 
modified for a Draine (1978J) interstellar radiation field. 
Hot CO production by HCO + recombinations is also in- 
cluded, as suggested by J. Black (2012, private commu- 
nication). 

For the X DR models, we used the grid by Meijerink| 
et al. ( 2007| ) that spans two different ranges in density 



(lTF^ 10 4 cm" 3 and 10 4 -10 6 5 cm -3 ) irradiation (F x 
1.6 — 160 ergs" 1 cm -2 ), and cloud size of 1 pc and 10 pc. 
We used a x 2 minimization algorithm to identify the best 
description of the 12 CO line fluxes measured with SPIRE 
and interferometric data (see Table IT]). 

As already mentioned in Section |3~1.1[ in addition to 
the CO lines detected in NGC 1068 byThe SPIRE FTS, 

13 Available at: |http: //www. strw.leidenuniv.nl/meije rink/grid/ 

14 Go = 1.6xl0 — ° erg cm"' s — 1 and the field is integrated from 
6 to 13.6 eV. 



with density ~ 10 cm -3 or ~ 10 b u cm -3 and incident 
X-ray flux in the range of 1.6-5.0 erg s" 1 cm~ 2 are able 
to reproduce the observed mid-J CO fluxes with a good 
X 2 value (see right panel of Figur e [7J . However, the 
RADEX models (see section 3.1.1) favor solutions with 
tih between 10 4 4 and 10 4 8 cm -3 as listed in Table [3j 

To further constrain the models, we have used the CO 
surface brightness for both PDR and XDR models, by 
summing over all of the modeled CO lines. For the CND, 
we have computed the total CO emission from the best 
fit XDR model presented in Figure [l0]and divided by the 
emitting area of the CND, assumed to be 4 " in diameter, 
and corrected by the filling factor of the RADEX best fit 
(see Table [3]) to derive the CO surface brightness of the 
CND. 

In Figures [8] and J9] we present the total CO sur- 
face brightness for botn PDR and XDR models, respec- 
tively, with the density region allowed by the RADEX 
models indicated by shaded regions. From Figure [51 
the predicted CO surface brightness in the CND by 
PDR models, which satisfies the density constrain of the 
RADEX model, is about two orders of magnitude below 
the value computed with the best fit model of the CND 
of 0.02 erg s -1 cm -2 sr _1 . On the other hand, Figure 
M shows that the predicted CO surface brightness in the 
CJND, for an XDR model with density of ~ 10 4 3 cnT " 



cm 



agrees well 
sr . 



and incident X-ray flux > 3 erg s 
with the best fit value of 0.02 erg s _1 cm 

We have also compared the observed parameters of the 
CND in NGC 1 068 with the XDR mod el results presented 
for Arp 220 in |Rangwala et al |p0TT| ). The XDR model 
of Arp 220 has been com puted using an update d version 
of the code described by Maloney et al.J ( 1996 ). The as- 
sumed physical conditions in Arp 220 coincide with those 
of NGC 1068: a hard X-ray (l-100keV) luminosity of L x 



Spinoglio et al. 




log n H [cm 



Fig. 7.— Reduced y 2 values of the models for the CND. Left: PDR models jWoIfire et al.||2010[ ) ; Right: XDR models jMeijerink" 
Spaans 2005 M eijerink et al.|2007| |. 




Fig. 8. — Total CO surface brightness (in units of 
erg s — 1 cm -2 sr _1 ) as a function of incident FUV flux and density 
for the PDR models. The left shaded region (in orange) shows the 
1<t range of the density of the best fit RADEX model of the SF-ring 
(extended region), while the right shaded region (in green) reports 
the area of the la range of the density of the best fit RADEX 
model of the CND (compact component). The red dashed line 
represents the measured surface density for the SF-ring, while the 
green dashed line represents the measured surface density for the 
CND of 0.02 erg s" 1 cm' 2 sr" 1 . 



~10 44 erg s _1 , a power law index of a =0.7 and an absorb- 
ing hydrogen column density of N# = 10 24 cm -2 , typical 
of Compton thick source s, like NGC 1068 . In th eir Figure 
7 (bottom, right panel) , [Rangwala et al. ( 2011[ ) show the 
predicted CO surface brightness as a function of density 
and radius of the emitting region. Our best fit value of 
0.02 erg s" 1 cm~ 2 sr -1 is in reasonable agreement with 



a density of ~10 cm" J and a radius of 100-200 pc. The 
radius that we associate to the CO emitting region in 
the CND is about 160 pc. We therefore conclude that 
the XDR origin of the CO CND emission in NGC 1068 
is indeed a plausible explanation. In Figure [10] we show 
the fit to the observed CO fluxes using the best fit XDR 
model with the adopted density and X-ray flux for the 
CND component. 




log n [cm ] 



Fig. 



9. — Total CO surface brightness (in units of 
1 ) as a function of incident FUV flux and density 
for the XDR models. The shaded region (in green) shows the lo~ 
range of the density of the best fit RADEX model of the CND 
(compact component). The dashed line (in green) represents the 
measured surface density. 

3.1.3. PDR and XDR models of the SF-ring component 

The low-J component of the SLED can be fitted either 
by PDR models with a FUV flux of G c > 10 2 and density 
of ~ 10 35 cm~ 3 , in rough agreement with the RADEX 
modeling results, or with models with lower FUV flux 
(G ~ 10 1 ) and a higher density of >10 4 cm' 3 , which are, 
however, in conflict with the RADEX results. Moreover 
(see the right of Figure 111, XDR models with a density 
~ 10 30 cm -3 and incident X-ray flux in the range of 0.5- 
2.0 ergs' 



cm 2 can fit the data. 



Sim ilarly to the case of the CND, described in Section 
|3.1.2[ for the SF-ring we have used the total CO emis- 
sion from the best fit PDR model and the area of 40" 
in diameter to derive the extended component surface 
brightness. For this latter case, we did not use any filling 
factor because we consider that the extended emission in 
the SF-ring is filling most of the beam. 

To further constrain the models, we use again the CO 
surface brightness. From Figure [8j the predicted CO 
surface brightness in the SF-ring, for a PDR model with 
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Fig. 10. — XDR model fitting for the compact comp onent, as- 
sociated with the CND. Models have been taken from [Meijerink| 
|fc Spaans| p005j and |Meijerink et al.| | |2007| l. The points follow 
the same color code as in f igurejbl The dotte d line shows the fit 
of the high-J CO lines with the XDR model of Hailey-Dunshcath 
|et al.| |2012), the dashed lines show the fit of the intermediate-J 
CO lines, with a XDR model with density n = 10 4 ' 25 cm~ 3 and 
X-ray flux of 9 erg s _1 cm -2 , once the high-J CO XDR model has 
been subtracted. The solid line shows the total modeled emission, 
which is compared to the observed data (see the text). 

density of ~ 10 3 cm~ 3 and incident FUV flux above Gq > 
10 2 , agrees well with the observed value of 2xl0 -5 erg 
s _1 cm" 2 sr _1 , reinforci ng this solution. We present the 
best fit model in Figure [12} 

3.1.4. Constraining the PDR and XDR models with 
luminosities 

In order to verify that the XDR models used to fit the 
compact (CND) and extended (SF-ring) molecular com- 
ponents in NGC 1068 are in agreement with the intrinsic 
X-ray flux emitted from the AGN, we have adopted an 
estimated intrinsic luminosit y in the range of Li mnbv 
~ 10 43 -10 43 5 erg s" 1 (see, e.g. [Colbert et al.|2002| ). Using 
the diameter of the CND of 4" and that of the SF-ring 
of 40", we derived an X-ray flux in the range of ~ 3-10 
erg s _1 cm" 2 and of ~ 0.03-0.11 erg s _1 cm -2 , respec- 
tively for the compact and extended components. 

The X-ray flux predicted for the XDR model of the 
CND component is in the range 1.6<Fjs: <10 erg s^'em" 2 
(see right panel of Figure [7]) , therefore the hypothesis of 
an XDR origin of the intermediate-J CO lines in the CND 
can be considered valid. The X-ray flux predicted for the 
XDR model of the SF-ring of F x = 0.5-2.0 erg s" 1 cm" 2 
is not in agreement with our estimates, being at least one 
order of magnitude higher, and therefore the XDR ori- 
gin of the lower- J CO lines in the extended region can be 
excluded, given the assumption in the above paragraph. 
Moreover, even if the X-ray flux would be available to 
excite the XDR, e.g. in the case of an underestimated 
intrinsic X-ray luminosity, an unknown part of it is ex- 
pected to be absorbed from the interstellar medium along 
the path, before reaching the SF-ring. The attenuation 
to the X-ray source as seen by the CO-emitting gas can- 
not in fact be quantified. Therefore the estimated flux 
could be considered as an upper limit. 

Similarly we estimated the far-UV (6-13.6 eV) flux 
from the AGN that illuminates the CND and the SF-ring. 
The adopted intrinsic far-UV luminosity of the AGN is 
~10 42 erg s _1 . This value was calculated using the in- 



trinsic far-UV continuum derived by Pier et al. (1994). 
At the distance of the SF-ring the far- U V flux is ^fuv = 
lO 5 G , whereas for the CND the F FUV = 10 2 4 G . 

This Fp\jv for the CND implies a high-density (riH > 
10 55 cm~ 3 ) component according to the PDR models 
(left panel of Figure ifh , which is not in agreement with 
the density estimate from the RADEX modeling. For 
the extended component, the UV flux from the AGN is 
lower than that predicted by the PDR models; however, 
in the star-forming ring it is likely that the UV emission 
from young stars also contributes to the interstellar UV 
radiation field. For comparison, th e far-UV flux from a 
young stellar cluster (age 5 Myr, see Spinoglio et al. 2005 ) 
of stellar mass 10 4 M is 10 2 G at 50 pc. This far-UV 
flux is an order of magnitude higher than that from the 
AGN and would be compatible with the predictions of 
the PDR models (see left panel of Figure 111. 

Taking into account all the above considerations, we 
therefore conclude that the most plausible explanation 
for the excitation of the CND is an XDR illuminated from 
the AGN, while that of the SF-ring is a PDR, mainly 
excited from the young stellar populations in the galactic 
arms. 

3.2. H2O emission 

We detected 6 H 2 transitions, 2 o-H 2 and 4 p-H 2 
lines, above a 3 a level within the SPIRE-FTS spec- 
trum. In Table [l] we list the detected lines, reporting 
their frequency, upper level energy E up /kb, and the mea- 
sured flux. We also listed the water lines detected with 
PACS (S. Hailey-Dunsheath 2012, private communica- 
tion). The line fluxes detected by PACS were measured 
with a Gaussian fitting of the line profile of the cen- 
tral spatial pixel (spaxel) and scaling using the standard 
point-source correction factors. The emission of all H 2 
lines detected in the SSW band of the SPIRE-FTS, as 
well as most of the PACS lines, are unresolved (diameter 
< 20", for the SPIRE lines, see Figure [TJ and < 10", for 
the PACS lines), indicating that water emission arises 
mainly from the CND region. Only one line, the p- 
H 2 2 n - 2 02 at 752.033 GHz was observed in the SLW 
band, with a beam size of ~ 33.4", and thus the emission 
of this transition may have an additional contribution 
from the starburst ring. However, our maps are not sen- 
sitive enough to disentangle the contribution of the star- 
burst ring from the CND. Moreover, the PACS map of 
the o-H 2 O2 12 -l i line at 179 ^m (1669.905 GHz) shows 
extended emission (S. Hailey-Dunshcath private commu- 
nication), most likely associated with the starburst ring. 
Thus, if there is some contribution from the extended 
emission in the central spaxel, the measured flux may 
overestimate the contribution of the nuclear flux. 

3.2.1. Excitation analysis 

In order to constrain the physical conditions of the wa- 
ter excitation, we built a grid of models using the Large 
Velocity Gradie nt (LVG) model in plan e parallel geome- 
try described in |Ceccarelli et al. ( 1998 ) varying the pa- 
rameters in the following range: n(H 2 )=10 3 -10 8 cm" 3 , 
T km = 10-1500 K, 7V(o-H 2 0) and A^(p-H 2 0) from 10 11 
to 10 17 cm" 2 . The molecular data were taken from the 



BASECOl]^] database ( |Dubernet et al.||2006| ) and we 

15 http://basecol.obspm.fr 
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Fig. 11. — Reduced \ 2 
the extended component, 
the model does not predict the lines of interest 



log n H [cm- 
values of the PDR jWolfire et al.|2010| and XDR models | |Meij 
Left: PDR models; Kigtit: ADK models. The line shaded regio: 



3.4 3.6 
log n H [cm 3 ] 



erink fc Spaans||2005[ |Meijerink et al.|2007| for 
ion ot the figure is where no model is available or 
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Fig. 12. — PDR/XDR model fitting for the extended component 
(SF-ring), together with the compact component, associated with 
the CND. Mode l s have been taken from Mcijcrink & Spaans j2005| ; 
|Meijerink et al.|l|2007|| . Notations as in b igure[b] Trie crosses (yel- 
low] represent the FTS measurements from the spectral map and 
the stars (white) s how the ground-based ob servations of the ex- 
tended component ||Kamenetzky et al,|201l| Wilson 2012, private 
communication). The thin solid lines show the fit of the compact 
compo nen t, which is the sum of the high- and mid- J models as in 
Figure [To| The dash-dotted lines show the fit of the extended com- 
ponent, with a PDR model with density n(H2) = 10 325 em~ 3 and 
log Go = 2.75, once the compact component has been subtracted. 
The thick solid line shows the total modeled emission, which is 
compared to the observed data (see the text). 



used the most recent collisional rate coefficients with H2 
Daniel et al.||2010l |2011h, calcu- 



(Dubernet et al. 2009 



lated lor temperatures between 5 K and 1500 K. The 
H 2 spectrum has been computed considering 45 levels 
for the ortho and para water (i.e., considering the levels 
up to excitation temperatures of ~ 2000 K) . The two wa- 
ter forms are treated as independent species. We used an 
ortho-to-para H2 ratio equal to 1 and an ortho-to-para 
H2O ratio of 1. The model includes the effects of the 
beam filling factor, so it computes the \ 2 f° r each col- 
umn density minimizing with respect to the source size, 
temperature, and density. Our model does not include 
the effects of radiative pumping. 
We fitted the o-H 2 and p-H 2 lines separately. We 



10.0 



a 

o 

• rH 

> 

u 

<D 

<si 

X> 

O 
"u 

o 



1.0 



0.1 



• N(p-H 2 0) = 8xl0 17 cm" 2 , n(H 2 ) = 5xl0 6 cm r 
A N(0-H 2 O) = 8xl0 17 cm" 2 , n(H 2 ) = 8xl0 7 cm" 3 



5 




f 



t 







100 



200 

E up /k B [K] 



300 



Fig. 13. — Ratio between the best fit model predictions and the 
observations as a function of the upper level energy. Triangles show 
the 0-H2O transitions and circles show the P-H2O transitions. The 
best fit models predict a column density of both ortho- and para- 
H2O of ~ 8 X 10 17 cm~ 2 and hydrogen desities in the range n(H2) 
= 5 x 10 6 - 8 x 10 7 cm" 3 . 



found that the observations of water lines are consis- 
tent with LVG model predictions for an emission size 
of 10" ± 2", a kinetic temperature of T kin ~ 40±10 K, a 
density n(H 2 ) from 4 x 10 6 to 8 x 10 7 cm , and a value 
of the N/Av varying between 5 x 10 14 and 10 16 cm~ 2 . 
The column density of each water form can then be 
obtained once we know the linewidth, which we set to 
100 km s _1 . This translates into a column density for 



both 0-H2O and P-H2O around 5x 10 
ure 
lines. 



16 



TO 15 



cm 



Fig- 



13 shows the best fit model for the 0-H2O and P-H2O 
The results of the LVG model predictions are sum- 
marized in Table|4] As can be seen in Fig. [T3j the best fit 
model predictions can reproduce most of the line fluxes 
reasonably well. Only for two 0-H2O lines the LVG model 
underestimates the observed line flux. However, the line 
flux of the 0-H2O 2i2-loi transition may be in fact over- 
estimated due to some extended emission that may come 
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from the starburst ring instead from the CND region (S. 
Hailey-Dunsheath priva te communication). As shown in 
previous sections and in Hailey-Dunsheath et al. (2012), 
the CO emission from the CND is associated with several 
components at different temperatures and densities. In 
the case of water we cannot rule out the presence of sev- 
eral components, but, due to the relatively low number 
of water lines, we did not attempt to fit the water lines 
with more than one component. 

Our results suggest a posteriori that the excitation of 
the water lines is not strongly affected by radiative pump- 
ing. Without the inclusion of radiative pumping, we can 
indeed fit all the detected lines within a factor of two. 
Moreover the water transitions in the SPIRE FTS spec- 
trum that should be more affected by radiative pump- 
ing, besides that of o-H 2 32i — 3i 2 , a re those of o-H 2 Q 
523 — 5i4 an d p-H 2 4 2 2 — 4 13 (see, e.g., |Gonzalez-Alfonso| 



et al.||2010[ ), which we do not detect, consistent with the 
LVG model predictions. On the other hand, we com 



pare our findings wit h the results of Mrk231 (Gonzalez- 
Alfonso et al. 2010), where all these three transitions 
have been detected and their fluxes can be reproduced 
only if radiative pumping is included. If we adopt the 
same line ratios of Mrk231 for NGC 1068, we would ex- 
pect to detect also these lines, originating from a higher 
level, at about a 10cr level (assuming the 3a sensitivity 
in the SSW FTS range of 10" 17 W m" 2 ). We conclude 
therefore that for NGC 1068 collisions seem to be the 
dominant excitation mechanism for water and that ra- 
diative pumping can be neglected. 

3.3. Molecular ions of H 2 + , OH+ and CH+ 

The molecular ions of H 2 + , OH + , and CH + have 
been detected in emission in NGC 1068. The 1115 GHz 
ground-state transition of H 2 + , as well as the 1033 GHz 
transition of OH + and that o f CH + , have also been de - 
tected in emission in Mrk 231 (Ivan der Werf et aTl[20To| . 
In contrast, these tra nsitions have all been detected in 
absorption in Arp 220 ( |Rangwala et al.|2011 ) and two of 
them (H 2 0+ at 11 15 GHz and OH+ afT0 33 GHz) in 



absorption in M82 (Kamenetzky et al 



2012) 



OH+, and 



To obtain the column densities of 
CH + , we assumed that the gas is in local thcrmody- 
namical equilibrium (LTE) and the emission is optically 
thin and originates in the CND region (diameter of 4"). 
Assuming that all levels are populated according to the 
same excitation te mperature, T ex , the column density is 
given by (see, e.g., Goldsmith fc Langer||1999[): 



N, 



mol 



47rFQ(T ex )exp(£ up /fcT cx ) 
hv A ul g u 



(1) 



where F is the integrated line brightness (in units of 
erg cm" 2 s" 1 sr" 1 ), v is the frequency of the transition, 
E up the energy of the upper level, Q(T ex ) is the parti- 
tion function at the excitation temperature T cx , and h 
and k are the Planck and Boltzmann constants, respec- 
tively. We approximated the partition function of CH + 
and OH + to Q ro t — kT ex /hB, where B is the rotational 
constant of the molecule. Einstein coefficients ^4 u i, up- 
per level degeneracy g u , and rotational constants have 
been obtained from the Cologne Database for Mo lecu- 
lar Spectroscopy (CDMS; puller et al.piOlj |2005[ ). For 
o-H 2 + the partition function has been obtained using 



the energy levels listed in the database. Both H 2 + and 
OH + present hyperfine structure, which is not resolved in 
the SPIRE FTS spectrum. Therefore, we used J2i M\9i 
to account for the several hyperfine components. We 
adopted an ortho-to-para ratio for H 2 + of 3. 

The estimated column densities are _/V(H 2 O + ) = (2.0 — 
5.9) x 10 13 cm" 2 , iV(OH+)=(1.0-1.2) x 10 13 cm" 2 , and 
A^(CH + )=(5.4— 7.5) x 10 13 cm -2 , for a range of excitation 
temperatures between 40 and 100 K. The uncertainty in 
the derived column densities is estimated to be around 
50%, which arises from the uncertainty in the flux cali- 
bration. The p-H 3 0+ J K = 3^ - 



2^ line at 364.7974 GHz 



has been de tected in NGC 1068 using the JCMT ( |Aalto 
et al. 2011). We recomputed the H 3 + column density 
using Eq. [TJ our adopted size of the CND region of 4", 
and the excitation temperature range of T ox ^ 40 — 100 K. 
We obtained iV(p-H 3 0+)~ (1.8 - 0.7) x 10 14 cm" 2 . As- 
suming an intermed iate ortho-to-para ratio of 1.5 



der Tak et al. 
(45" 



2008 

x 10 ; 



van 



the total column density of H 3 + 



IS — (4.t> - 1. 

Considering the sp ectrum of Arp 220 presented in 
Rangwala et al. (2011 ), we notice that in NGC 1068, the 
ratio of JN (H 3 O^77N(H 2 O + )~10, whereas this ratio is 
reversed in Arp 220. Similarly, the inferred column den- 
sities indicate that N(H 3 0+)/N(OH+) > 1 in NGC 1068, 
while in Arp 220 the OH + column density is at least an 
order of magnitude larger than that of H 3 + . So the 
relative abundances of the molecular ions appear to be 
very different in the two galaxies, suggesting some sub- 
stantially different conditions. 

We compared o ur estimated column den sities with the 
predictions of the |Hollenbach et al. (2012) PDR models 
for different FUV incident fluxes and cosmic ray ioniza- 
tion rates. For the FUV radiation field produ ced by the 
AGN in the CND region (see Section 3.1.4) the order 
of magnitude of the observed column densities of H2 + 
and OH + are comp atible with a low density medium (see 
Fig. 9 and Fig. 10 of|Hollenbach et al.| (f2012l)) . According 



to their models, the column density of M 3 + is approx- 
imately constant with a value of 10 14 cm" 2 , in agree- 
ment with our estimated value. Combining the observed 
value of the N(OH+)/N(H 2 0+) m tio and the value of 
A^(OH + ), we estimated a range of values of logio[£/^] 
between —15.5 and —16.5, where £ is the cosmic ray ion- 
ization rate and n is the hydrogen nucleus density, and 
we derived A^h > 2 x 1 20 cm -2 (see Figs. 10 to 13 of 



Hollenbach et al.| ( |2012) ). Therefore PDR models which 
include cosmic rays are in agreement with the observed 
column densities of these molecular ions. 

3.4. Hydrogen fluoride 

Herschel spectroscopic observations have detected for 
the first time the J = 1-0 transition of hydrogen fluo- 
ride (HF) at 1232 GHz in the local universe, and re- 
vealed the ubiquitous nature of this molecule in the in- 
terstellar medium (ISM) of the Milky Way and of a few 
extragala ctic sources, such as the local AG N/ULIRGs 
Mrk 231 (Ivan der Werf et~aT1[2010]), Arp 220 ( |Rangwala 
et al.||201T| ), M82 ( |KamenetzlcFetal~||2012|) a nd in the 
(Joverleat quasar at z=2.56 ( |Monje et al.||201l[ ). 

This transition is generally observed in absorption, as 
expected, due to its very large Einstein coefficient, Aio = 
2.42 x 10" 2 s _1 . Only an extremely dense region, with a 
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strong radiation field, could generat e enough excitation 



to yi eld an HF feature in emission (Neufeld & Wolfire 
20091. In NGC 1068 it has been detected in emission (at 



a S/N > 10), similarly to Mrk 231 and in contrast with 
Arp 220, where it has been detected in absorption. 

We estimated the column density of hydrogen fluoride 
using Equation [T] assuming that the HF emission origi- 
nated in the CND region (diameter ~ 4"). We adopted 
the kinetic temperature Tki n — 100 K derived for the CO 
gas associated with the CND region (see Sec. 3. 2.1) as an 
upper limit to the excitation temperature. The values of 
the upper level degeneracy, g up = 3, and the rotational 
constant (_B=616365 MHz) to calculate the partition 
function at a given temperature, Q TOt — kT cx /hB, have 
been ta ken from the Jet P ropulsion Laboratory (JPL) 
catalog ( |Pickett et al]|1998| . For T ox = 40 — 100 K, we 
found a column density of HF of iV(HF)~ 7 x 10 12 cm' 2 . 
This value is one order of magnitude lower than the 
N(RF) value found in M82 (Kamenetzky et al. 2012), 
where hydrogen fluoride is seen in absorption. Calcula- 
tions with RADEX indicate that the optically thin ap- 
proximation holds for a wide range of physical conditions, 
with t<1 for temperatures between Tkin = 20 — 100 K, 
densities n(H 2 )~ 10 3 - 10 6 cm 3 and column densities 
of 10 n -10 14 cm' 2 . 

3.5. Atomic lines 

In this section, we analyze the origin of the three 
atomic fine-structure lines detected with the SPIRE 
FTS, namely, the two lines from neutral carbon and the 
[Nil] 205 fim line. First we atte mpt to use the composit e 
photoionization model used by Spinoglio et al. (2005), 
which successfully fitted the overall U V to tar-lR spec- 
trum of NGC 1068 reproducing the line fluxes within a 
factor 2 on average. This model was composed of an 
AGN component and a starburst component, which in- 
cluded contribution from PDR clouds, as the integration 
was allowed to run until the gas temperature in the cloud 
cooled down to T=50 K. We compare in the next section 
the predictions of the atomic lines detected by the SPIRE 
FTS using the s ame p hotoionization model reported in 
Spinoglio et al. ( 2005 ) , which were not reported in that 
work. 

3.5.1. [Nil] and [CI] emisson 

The observed [Nil] 205 /im line is roughly consistent 
with the p hotoionization model presented in | Spinoglio] 
et al.| ( |2005[ ). The predicted ratio [Nil] 205 /xm/122 //m is 
~ 0~47 wEile the observed value is 18.8/30.5= 0.62±0.09, 
using the ISO-LWS ob servation of the [Nil] 122 //m line 
(Spinoglio et al. 2005). Taking into account that the 
two fines have been observed by two different instru- 
ments, with different apertures, this result is satisfac- 
tory. We can therefore conclude that the diffuse ionized 
emission traced by the [Nil] ion is in agreement with be- 
ing excited by the starburst component in the SF-ring 
of NGC 1068. The observed [CI]369 /zm/[CI]609 /mi line 
ratio is 1.65±0.30, while the ratio predicted by the pho- 
toionization model is about 3.8, indicating that the emis- 
sion of the line at 609 fim could be more extended than 
the FTS beam. 

On the other hand, the intensity of the [CI] lines pre- 
dicted by the photoionization model are weaker by more 



than an order of magnitude compared to the observed 
values r elatively to the [ CII]15 8 emission, as mea- 
sured in |Spinoglio et al." (2005). The predicted ratio of 
[CI] 369 /zm/[CIIJ158 fim is in fact - 10" 3 and the pre- 
dicted [CI]609/[CII]158 fjtm fim is - 3 x 10~ 4 , while the 
observed ratios are 9.7 x 10 -3 and 5.9 x 10 -3 , respec- 
tively. This discrepancy between data and models could 
be due to the poor ability of current models to reproduce 
the [C I] emission line intensities (see, e.g., |Rollig et al.| 
2007). 



3.5.2. Using Neutral Carbon emission to estimate the gas 
temperature 

We derive here the kinetic temperature of the molec- 
ular gas from the ratio of the intensities of the neutral 
Carbon lines ([CI] at 369 and 609 /im), when local ther- 
modynamic equilibrium (LTE) is assumed. The Boltz- 
mann equation can be written as 



T k = 



-E 



21 



gin 2 

92TII 



(2) 



where 7\. is the kinetic temperature of the gas, E21 = 
E2 —Ei, n\ and tj 2 are the level populations and g\ 
and </ 2 their statistical weights. The integrated flux in 
a line is simply cx A^n^^. Therefore the kinetic 
temperature can be written as 



T, 



-E21 



1 1 9i 

In — 



-Aio^io-Fbi 



92 Ai\V2\F_ 



10 



(3) 



By substituting the observed values of the fluxes of the 
[CI] lines at 492.16 and 609.34 GHz of Tableland propa- 
gating the uncertainties, we obtain a kinetic temperature 
of the gas of T^ = 22.5 ± 2.4 K. The kinetic temperature 
derived from the two [CI] lines is much lower than for ei- 
ther the CND or extended starburst ring components, 
which again argues that the CI emission is much more 
extended than either of these. 

3.6. Comparing the derived masses with previous work 

We discuss here our results in the context of the vari- 
ous molecular observations of NGC 1068 in the literature. 
In particular our radiation transfer models indicate the 
presence of two major components responsible for the CO 
(J«p _ 13) excitation: the first one is a compact compo- 
nent, with a diameter of ~ 4" (300 pc), associated with 
the CND with a density of n(H 2 ) ~ 4xl0 4 cm -3 , kinetic 
temperature of T-90 K and mass M(H 2 ) ~ 2.4x 10 7 M , 
and the second is an extended component, with a diam- 
eter of ~ 40" (3 kpc), associated with the SF-ring, with 
a density of n(H 2 ) ~ 7xl0 2 cm' 3 , kinetic temperature 
o f T~116 K and mass M (H 2 ) ~ 3.5x10 s M w . 



Schinnerer et al. (2000) have used interferometric CO 
(1-0) and (2-1) observations to derive a mass of the CND 
("ring" in CO (2-1), in their terminology, of ^200 pc 
in radius) of M(H 2 ) ~ 5xl0 7 M©, while they measure 
a mass of the "spiral arms" (our SF-ring) of M(H 2 ) ~ 
6.8 xlO 8 M Q . These two mass estimates are a factor 2 
times larger than our estimates; however, this discrep- 
ancy is fully justified by the difference in the gas excita- 
tion, as the CO (1-0) and (2-1) lines are mapping lower 
temperature gas. 
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We then compare our findings with the temperature of 
T>145 K and mass derived from the H 2 pure rotational 
Rigopoulou et al7| (|2002[ ); from the detec- 



transitions of 
tion of the S( 
of M(H 2 ) ~ 10 
x 27" ( 



line (Lutz et al 

M, 



a warm gas mass 



1R 



is derived in the SWS beam of 14" 



de Graauw et al. 1996). This is consistent with 



our findings, given the difference in beam size between 
the two spectrometers. 

3.7. Observability of submillimetre lines at high redshift 

Using NGC 1068 as a template to predict the submil- 
limetre spectrum of higher redshift galaxies, we have es- 
timated the expected line fluxes as a function of redshift, 
rescaling the total infrared luminosity of NGC 1068 to the 
value of L/^=10 12 Lq. We show in Figure [l4[a) the ob- 
served fluxes of atomic and molec ular lines of NGC 106 8 
from the ISO-LWS spectrometer (Spinog lio et al.||2005|), 
from the PACS spectrometer (Hailcy-Dunsheath et al.| 



2012 1, and those from this work, and compare them with 
the predicted fluxes at redshift of z=0.1, 0.2, 0.3, 0.5 and 
1, 2, 3, 5. We also show in the figure the 5cr 1 hour sensi- 
tivities of future observing facilities, such as CCAT and 
SPICA-SAFARI and compare them with the sensitivi- 
ties o f the H erschel spectrometers. We refer to Spinoglio 



et al. ( 2012J) and references therein for a brief description 
ot SPltJA and CCAT and for the details of the expected 
sensitivities of the foreseen spectrometers at their focal 
planes. We also show in Figure [l4]^b), for comparison 
with NGC 1068, the observed and predicted submillime- 
tre spectra of other important local template galaxies, 
the prototypical starburst galaxy M82 and the ULIRG 
Arp220, t hat have been observed with the SPIRE spec 



trometer ( |Kamenetzky et al.|2012 Rangwala et al.|2011 
respectively). 



respectively) 

It is clear from Figure 14 that the intermediate to high- 
J CO lines, as well as the two [CI] lines, are weaker by 
1-2 orders of magnitude compared to the brightest far-IR 
fine structure lines of [CII] and [Nil]. However, their pow- 
erful diagnostic potential, in terms of detecting through 
XDR regions the effect of an AGN in the host galaxy, 
combined with their long wavelengths, makes these lines 
very attractive for high redshift spectroscopic cosmologi- 
cal surveys from ground-based telescopes. The expected 
sensitivity of future spectrometers at the focal plane of 
large submillimetre telescopes, such as CCAT, will be 
able to detect an object like NGC 1068, M82, or Arp220 
with a luminosity of L=1O 12 L at a redshift of z=0.2-0.5. 

4. SUMMARY AND CONCLUSIONS 

We summarize here the results of this work. The first 
complete submillimetre (190-670 /im) spectrum of the 
Seyfert type 2 galaxy NGC 1068 reveals the full sequence 
of CO pure rotational lines from J up =4 to J up =13. The 
radiation transfer analysis of these lines shows the pres- 
ence of two physically distinct components: the first one 
originating from the circum-nuclear disk (CND) of few 
arcseconds in diameter (~ 4") and the second one ex- 
cited in the star forming ring (SF-ring) with a diameter 
ten times larger (~ 40"). These results indicate a ki- 
netic temperature of CO of lkj n =100 K and 127 K, a gas 
density of ri(H 2 )=10 45 and 10 29 cm -3 and a derived 
molecular hydrogen mass of M(H2) ~ 2.4xl0 7 M and 
M(H 2 ) ~ 3.5x10 s M Q , for the compact and extended 
regions, respectively. 



The comparison of the observed CO line intensities 
with predictions of photodissociation (PDR) and X-ray 
dominated regions (XDR) models shows that the circum- 
nuclear disk emission can be modeled equally well by 
both types of models, while the CO lines in the star- 
forming ring can be modeled by a photodissociation re- 
gion only. However some observational constraints, such 
as the total CO surface brightness and the required radi- 
ation field, indicate that the most plausible explanation 
for the CO excitation of the CND is an XDR with density 
of n(H 2 ) ~ 10 4 cm -3 and X-ray flux of 9 erg s _1 cm~ 2 , 
consistent with the AGN illumination. In contrast, the 
excitation of the SF-ring component is due to PDR emis- 
sion originating from the young stars/HII regions in the 
spiral arms. 

The water lines that we have detected with SPIRE, 
together with those observed by PACS (S. Hailey- 
Dunsheath 2012, private comm.), have been modeled 
with an LVG model to constrain the physical conditions 
of the water excitation. We have found that the kinetic 
temperature is Tfci„=40 K, the molecular hydrogen den- 
sity is n(£f 2 ) = 4 x 10 6 — 8 x 10 7 and the column density 
is of order N(H 2 0)=8 xlO 17 cm~ 2 for both water forms. 

The computed column densities of the molecular ions 
detected (H 2 + and OH + ) are in agreement with PDR 
models that include cosmic ray ionization. 

The fundamental rotational transition of HF has been 
detected in emission in NGC 1068 and we infer a column 
density of N(HF)~ 7.5 x 10 12 cm" 2 . 

For the two [CI] transitions, we derived a kinetic tem- 
perature of 22.5±2.4 K in LTE approximation, which 
is much lower than the temperatures traced from the 
intermediate-J CO molecular gas. 

The molecular masses that we derived from our analy- 
sis are in good agreement with the masses estimated from 
both CO interferometric measurements of low- J lines and 
mid-infrared H 2 emission lines. 

Finally we show that the intermediate-J CO and [CI] 
lines in galaxies with L~ 10 12 L Q can be observed from 
planned and future ground-based and space telescopes up 
to redshift of z ~ 0.5, making their diagnostic power an 
important tool to study galaxy evolution at intermediate 
redshift. 
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Fig. 14. — Left panel: Comparison between the observed far-IR/submillimetre lines of NGC 1068 with the predictions of the same lines 
as a function of redshift of a similar galaxy with an intrinsic luminosity of L=10 12 Lq. The stcllar/HII region lines are shown as filled 
triangles (in red), the low-ionization/PDR lines are shown as filled squares (in green), while the molecular lines are shown as open circles 
(in magent a). For comparison, we overplot the 5cr (in 1 hr) sensitivity threshold of Herschel-PACS, Herschel-SPIKE, SPICA-SAFARI and 
CCAT. See |Spinoglio et al.| ( [2012] ) for more details on the assume d sensitivities of the futur e SAFARI and CCAT ins truments. Right vane l: 
Same diagram tor other two templates: the starburst galaxy M82 l |Kamenetzky et al.|2012] | and the ULIRG Arp220 (Rangwala ct al. |2011| |. 
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TABLE 1 

Fluxes of the molecular and atomic lines detected in the SPIRE FTS spectrum of NGC 1068 and compared to 

LITERATURE DATA 

species line id. rest. freq. E up /Kg Beam Flux la Stat. Uncert. lcr Tot. Uncert. Ref. 

FWHM 



(GHz) (K) (") (10~ 17 W m" 2 ) 



compact source 





J = 1-0 




115.271 


5.53 


1 x 0.8 


.050 


.001 


0.01 


1 




J = 2-1 




230.538 


16.60 


1 x 0.8 


.360 


.003 


0.07 


1 




J = 3-2 




345.796 


33.19 


1 x 0.8 


2.36 


0.34 


0.50 


1 




T O 

j — y-o 




1 nQfi 01 


O/l Q QQ 

z4o.oo 


lo. I 


1 K A 1 

10.41 


n 1 
U.10 


1.00 







T — in n 

j — J- u- y 




11^1 QR^ 




17 1 


14 32 a 




1 


2 




J = 11-10 




1267.014 


364.97 


17.6 


12.80 


0.15 


1.29 


2 




J = 12-11 




1381.995 


431.29 


16.9 


10.97 


0.18 


1.11 


2 




J = 13-12 




1496.923 


503.13 


16.8 


6.11 


0.27 


0.67 


2 


o-H 2 


3l2 — 3()3 




1097.365 


215.1 


18.7 


4.15 


0.03 


0.42 


2 




321 — 3l2 




1162.912 


271.0 


17.1 


5.80 


0.14 


0.60 


2 




221 — 2i2 




1661.008 


159.8 


12.7 


4.8 




1.6 


3 




2l2 — 1()1 




1669.905 


80.1 


12.7 


14.3 




4.3 


3 




303 — 2 12 




1716.769 


162.5 


12.3 


6.9 




2.2 


3 




4l4 — 3()3 




2640.474 


289.2 


8.0 


5.1 b 




1.6 


3 




221 — llO 




2773.978 


159.8 


7.6 


6.5 




2.1 


3 


P-H2O 


2n — 2o2 




752.033 


136.9 


33.4 


3.42 


0.07 


0.35 


2 




2()2 — 111 




987.927 


100.8 


19.8 


6.09 


0.11 


0.62 


2 




In — Ooo 




1113.343 


53.4 


18.7 


3.49 


0.05 


0.35 


2 




220 —2i.i 




1228.789 


195.9 


17.6 


5.30 


0.09 


0.54 


2 




322 — 3l3 




1919.359 


296.8 


9.8 


2.4 




0.8 


3 




3l3 — 2 02 




2164.132 


204.7 


11.5 


6.9 




2.2 


3 


o-H 2 CT 


h In - Ooo J= 


3/2-1/2 


1115.204 


53.5 


18.7 


2.13 


0.05 


0.22 


2 


o-H 2 CT 


h In - Ooo J= 


1/2-1/2 


1139.561 


53.5 


17.1 


1.98 


0.09 


0.22 


2 


OH+ 


112 — 012 




1033.118 


49.6 


18.7 


4.34 


0.16 


0.46 


2 


CH+ 


J = 1-0 




835.079 


39.6 


35.0 


3.16 


0.08 


0.33 


2 


HF 


J = 1-0 




1232.476 


59.1 


17.6 


2.16 


0.09 


0.23 


2 


[Nil] 


3 Pi - 3 Po 




1461.134 


6.77 


16.8 


54.97 


1.42 


5.67 


2 


extended source 


12 CO 


J = 2-1 




230.538 


16.60 


31.0 


6.41 


.015 


1.3 


4 




J = 3-2 




345.796 


33.19 


43.4 


19.4 




3.4 


5 




J = 4-3 




461.041 


55.32 


43.4 


21.69 


0.66 


2.73 


6 




J = 5-4 




576.268 


82.97 


43.4 


26.80 


0.43 


3.30 


6 




J = 6-5 




691.473 


116.16 


43.4 


24.73 


0.24 


3.02 


6 




J = 7-6 




806.652 


154.87 


43.4 


21.06 


1.30 


2.88 


6 




J = 8-7 




921.800 


199.11 


43.4 


19.78 


0.64 


2.50 


6 


[CI] 


3 Pi - 3 Po 




492.161 


23.62 


43.4 


12.81 


0.13 


1.57 


6 


3 P 2 - 3 Pi 




809.342 


62.46 


43.4 


21.08 


1.20 


2.84 


6 


[Nil] 


3 Pi - 3 Po 




1461.134 


6.77 


43.4 


187.9 


8.2 


24.3 


6 



References. — (1): Interferometric observations from [Krips et al.| l |20iT[ | , adding together the emission from the E-knot and 
the W-knot and using the conversion: /(W m~ 2 ) = 3.33 Xl0~™ i/gh» Sco ( Jy km s" 1 ); (2): This work SPIRE-FTS deep pointed 
spectrum; (3): S. Hailcy-Dunshcath, 2012, private communication; (4): From CSO Z-Spec observations of Kamcnctzky et al.| l |201l[ , 
using the same conversion as for (1); (5): C. Wilson, 2012, private communication; (6): This work, SPIRE-FTS line maps 

Note. — a This line is possibly contaminated by the 0-H2O 3i2 — 221 at 1153.127 GHz. However, the line fitting results indicate 
that the CO line is dominant; therefore we assign the whole flux to the CO line; This line is blended with t he CO (23—22) transition, 
so to obtain the flux of water we subtracted the flux of CO (23—22) predicted by the model presented by |Hailey-Dun shcath ct al. 
120121. 
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TABLE 2 

RADEX Model Parameters and Ranges 



Parameter 


Range 


# of Points 


Tkin [K] 


10 1.7_ 10 3.5 


90 


n(H 2 ) [cm- 3 ] 


10 10 - 10 5 


90 


A^CO [cm" 2 ] 


10 16 _ 10 20 


90 


AV [km s" 1 ] 


1.0 


fixed 



Note. — All parameters are sampled evenly 
in log space. 



TABLE 3 

Multi-Component Radiative Transfer Results 



Parameter Median lcr Range ID Max 4D Max 



compact source 



T kin [K] 


92 


74 - 


120 


87 


100 


n(H 2 ) [cm- 3 ] 


1Q 4.6 


io 4 - 4 - 


1Q 4.8 


1 4.6 


1 4.5 


N co [cm- 2 ] 


10 19.5 


10 19.3 _ 


1 19.7 


1Q19.5 


1Q19.6 




10 -0.9 


10-1.0- 


1Q-0.9 


1Q-0.9 


1Q-0.9 


P [K cm- 2 ] 


10 6.6 


106-5 _ 


1Q6.8 


1Q6.6 


1Q6.6 


{Nco) [cm- 2 ] 


10 18.6 


io 18 - 4 - 


1 18.7 


10 18.6 


1Q18.6 


M(H 2 ) [M ] 


io 7 - 4 


io 7 - 2 - 


1 7.5 


io 7 - 4 


io 7 - 4 



extended source 



T kin [K] 


116 


80 


- 166 


121 


127 


n(H 2 ) [cm" 3 ] 


1Q2.8 


1 2.3 


- io 3 - 3 


1 2.9 


1 2.9 


N co [cm" 2 ] 


1Q18.8 


1 18.3 


- ioi 9 - 3 


1 18.7 


1 18.7 


*A 


io-i- 2 


10-1.3 


- IO" 10 


io-i- 2 


io-i-i 


P [K cm" 2 ] 


io 4 - 7 


1Q4.3 


- io 5 - 2 


1 4.4 


io 4 - 4 


(Nco) [cm- 2 ] 


1Q17.7 


1Q17.2 


- ioi 8 - 2 


1 17.6 


io 17 - 6 


M(H 2 ) [M ] 


10 8.5 


1 8.0 


- 10 90 


1 8.5 


10 8.S 



Note. — ID max refers to the maximum value of the integrated pa- 
rameter distribution. 4D max refers to the value of that parameter at 
the best fit solution. 



TABLE 4 

The best fit model parameters of 
H 2 Q 



Parameter Value 

Temperature 40 ± 10 K 

density 4 X 10 6 - 8 X 10 7 cm" 3 

diameter 10" ± 2" 

W(o-H 2 0) 2 x IO 47 - 10 18 cm" 2 

JV(p-H 2 0) 5 x lQi 6 - IO 28 cm- 2 



